The bosonic birthday paradox 
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We motivate and prove a version of the birthday paradox for k identical bosons in n possible modes. If the 
bosons are in the uniform mixed state, also called the maximally mixed quantum state, then we need k ~ \fn 
bosons to expect two in the same state, which is smaller by a factor of v2 than in the case of distinguishable 
objects (boltzmannons). While the core result is elementary, we generalize the hypothesis and strengthen the 
conclusion in several ways. One side result is that boltzmannons with a randomly chosen multinomial distribu- 
tion have the same birthday statistics as bosons. This last result is interesting as a quantum proof of a classical 
probability theorem; we also give a classical proof. 
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The traditional birthday paradox says that given a calen- 
dar with n days, there is a significant chance (bounded away 
from 0) that a room with £l(^/n) people with uniformly ran- 
dom birthdays has two with the same birthday. Aaronson and 
Arkhipov 11J] discuss the same paradox for randomly chosen 
bosons. Here we present a different treatment of the same 
problem. In fact we will present two "paradoxes". The first re- 
sult (which Aaronson and Arkhipov derived, in a less general 
form) is that although bosons prefer to have the same birth- 
day, they have the same asymptotic behavior in the birthday 
problem, up to constant factors, as distinguishable particles 
(boltzmannons). The second result is that they have exactly 
the same behavior, non-asymptotically, as n i.i.d. boltzman- 
nons whose common distribution is a randomly chosen point 
in the simplex of all distributions on n configurations. This 
leads to an interesting result in classical probability with a 
quantum probability proof. 

We assume that the Hilbert space for one particle is 3tf = 
C". We assume a self-adjoint birthday operator 

B : ,ye ^r,ye 

with eigenvalues 1,2, ... ,n in some basis. The Hilbert space 
of k bosons is then the symmetric power 

s*(jr)sc((*)), 

using the multiset coefficient notation 



def (n + k — 1 

k 



(1) 



In the terminology used for identical particles, the states of a 
basis of J4? are called modes. 

In the traditional version of the classical birthday problem, 
we assume the uniform distribution on all n choices of the 
birthdays of the k people. The uniform distribution /ViifPO 
on any finite set X can be characterized in either of two ways: 
It is the unique distribution with the most entropy, log \X\; and 
the unique distribution with the most symmetry, Sym(X). 
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We will consider an analogue of the uniform distribution 
for a quantum system with a Hilbert space the mixed 
state p un if (Jf?) whose density matrix is the scaled identity on 
Jff. Like the classical state /i un if, the quantum state p un ;f (Jf?) 
is the unique state on with the most entropy, log dim 
and the unique state with the most symmetry, U (Jf?). More- 
over, Punif(^) is the unique state that yields the distribution 
Humf(X) for any complete measurement that takes values in a 
setX. 

We will use the uniform state p un ;f = p U mf{S k {J4f)) on the 
joint Hilbert space of k bosons. Then, the measurement S (B) 
of all birthdays of p un ;f yields the uniform distribution /x un ;f 
on configurations of k unlabelled people with n possible birth- 
days. (It is also standard to refer to unlabelled balls in labelled 
boxes, but we will stick to the birthday metaphor.) Moreover, 
this particular uniform state can be justified using less sym- 
metry than the largest available unitary group U(S k (J^)): 

Proposition 1. The state p un if on S k (J4?) is the unique state 
which is invariant under the unitary group U{J^C). 

Proof. Suppose that p is a U{M' )-invariant state on S (J%? ), 
i.e., a t/(j£°)-invariant density operator. Schur's Lemma says 
that if V is an irreducible complex representation of a group 
G, then every G-invariant operator on V is proportional to the 
identity. Thus it is sufficient (and also necessary, if either V 
is unitary or G is compact) for V to be irreducible. It is a 
standard fact of representation theory §6.1] that S k (C" ) is 
an irreducible representation of GL(n,C). It is another stan- 
dard fact §26.1] that GL(«,C) and U(n) have the same 
irreducible representations, since the former is the complexi- 
fication of the latter. □ 

This symmetry implies that p un ;f is the (7(J^)-average of 
any state, since such an average must be invariant with respect 
to the action of U(Jf). 

Corollary 2. Putting k bosons in any state <T on S k (Jif), and 
then applying a Haar-random unitary matrix in U (J^f) yields 
the state p un ;f. 

Aaronson and Arkhipov consider such an average for a par- 
ticular choice of <7, where a is the pure state 

|y> = |l,2,3,-,*) 
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in which the k bosons are in distinct modes (which requires 
k<ri). Another choice considered below is 

|vr) = |l,l,l,...,l> 

in which the bosons are all in the same mode. There are many 
choices for a, but Corollary |2] says that they all become the 
same when they are averaged. 

We will now look at the asymptotics of y'-fold birthdays 
in Punif- We will use the notation f(n) ~ g(n) to mean that 
f(n)/g(n) — > 1, or equivalently that/(«) = g(n)(\ +o(l)). 

Theorem 3. Suppose that there are k bosons with n modes, 
suppose that they are in the uniform state p un if, and suppose 
that k ~ cn^~^'j as n — > °° for some integer j>2 and some 
constant c > 0. Then the number of j-fold birthdays converges 
in distribution to a Poisson random variable with mean c 1 , 
while the number of (j + \ )-fold-or-more birthdays converges 
toO. 

This is the same asymptotic answer as in the case of boltz- 
mannons, except that the mean in that case is c 7 / j\. In fact, 
our argument in the case of bosons is very similar to a standard 
argument in the case of boltzmannons. 

Proof. Recall that the joint measurement S k {B) of all of the 
birthdays yields the uniform distribution on k unlabelled peo- 
ple among n calendar days. The probability that the first birth- 
day has at least j + l people is 



n 

k-j-l 



{n + k)i+ l 



for fixed j and n,k~^> I. Taking k = 0(n^ '■"■') and summing 
over all n days, the expected number of (j + l)-fold-or-more 
birthdays is 0(n -1 '-'), which vanishes as n — > °°. 

Meanwhile the probability that the first I days each have at 
least j people is 



n 

k-ji 



n 



k>< 



{n+k)-i 11 



where the approximation holds for fixed j and I and n,k^> 1. 

Summing over all (") ~ |y choices of the £ days, we obtain 
that if X is a random variable representing the number of j- 
fold birthdays, then 



X 



C J1 



So in the limit, the £th factorial moment is c^, which the same 
answer in the limit as a Poisson random variable with mean 
cK To conclude the argument, the Poisson distribution is de- 



termined by its moments. 



□ 



The calculation for the narrow question of the probability 
of at least one repeated birthday is simpler. The probability 
that all of the birthdays are distinct is 



k-l 

n 



,-*7» 



as long as k = o(« 3//4 ). The approximation is established by 
taking the logarithm of both sides and then applying the Taylor 
series estimate 
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-2x + 0{x 3 ). 



Corollary 4. For n modes, we need k ~ %/«ln2 bosons to 
expect a repeated birthday with majority probability. 



This differs by only a constant factor from the k ~ ^2«ln2 
people needed to expect a repeated birthday in the classical 
birthday problem with distinguishable people. 

Remark. We should say something about independent but 
non-uniform bosons. The notion of independence for bosons 
is subtle. One reasonable and widely used notion is to first 
choose a distribution jj. for the birthdays of one boson, and 
to model it by a diagonal density matrix in the birthday ba- 
sis. Then there is a unique distribution on k bosons such 
that if A: — 1 of the bosons are fixed, the conditional distri- 
bution of the last one is given by jj.. This distribution is also a 
thermal state, also known as a Maxwell-Gibbs state, for non- 
interacting bosons. It was discovered by Bose and Einstein 
that under fairly mild assumptions on fl, almost all of the 
bosons have the most likely birthday. This paradox is com- 
monly known as Bose-Einstein condensation. 

Corollary [2] implies an interesting second model for the 
joint distribution of birthdays of k bosons. 

Theorem 5. The joint birthday distribution ofk bosons in the 
uniform state p un if is identical to the average ofk i.i.d. boltz- 
mannons, if their common distribution is given by a uniformly 
random point in the simplex of distributions on the n birth- 
days. 

By combining with the induced uniform distribution on the 
birthday measurement, we obtain a corollary of Theorem [5] 
that equates two distributions in classical probability. 

Corollary 6. Consider a town in which all families first agree 
to have children according to a common distribution on the 
days of the year, which itself is chosen uniformly from the sim- 
plex of all distributions. Then the children's birthdays behave 
as if the children were unlabelled, i.e., if we make a table that 
only gives the number of children born on each day, then all 
such tables are equally likely. 

In other words, the uniform average of all multinomial dis- 
tributions on multisubsets of size A: in a set of size n, is the 
uniform distribution on multisubsets. 

Proof of Theorem\5\ Recall that the Hilbert space of k boltz- 
mannons is y?® k . Consider the state a = (| y) {y\)® k , first for 
some fixed choice of | £ Jff. This a yields independently 
distributed birthdays for the k boltzmannons, and the distribu- 
tion of each one is given by the measurement of one copy of 
| y) . Meanwhile, a is evidently a pure symmetric state, which 
means that these boltzmannons are also bosons. By Corol- 
lary |2j the average of all choices of <T, with respect to Haar 
measure on U(Jt?), is the bosonic state p un if. 
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The Haar distribution of 1 1//), or equivalently one column of 
a matrix in U(n), is given by Haar measure on the manifold 
of pure states CP"~ 1 . The induced distribution of the birthday 
measurement is given by the moment map 

m: CP"' 1 ->A„_i 

to the simplex of distributions on n configurations UA §6.4]. 
This moment map preserves normalized measure |]2l §6.6]. 
Thus a random choice of a amounts to a random distribu- 
tion on each birthday, drawn uniformly from the simplex of 
distributions. This establishes the claim of the theorem. □ 

Theorem[5]yields a quantum proof of a classical probability 
result, Corollary [6] We also obtained a classical proof of the 
same result. 

Classical proof of Corollary® The argument uses a variation 
of the stars-and-bars notation for multisets [0] that is also used 
to prove the identity ([T). Namely, we write a star for each 
of the k children, with n — 1 separating bars between the n 
calendar days. For example, if there are k = 4 children and n = 
6 birthdays, then one possible choice for all of the birthdays is 

★★1*11*11, 

in which two children are born on the first day, one on the sec- 
ond day, one on the fourth day, and none on the other days. 
We first choose locations of n— 1 bars independently and uni- 
formly on the unit interval / = [0,1]. This separates the inter- 
val into n subintervals of length 

Pl+P2 + --- + Pn = h 

and we claim that the lengths of these subintervals are given 
by a uniformly random point in the simplex of distributions. 



(Because, if we first take the bars to be numbered, they are 
distributed according to uniform measure on [0,1] . Then, 
erasing the numbers yields the quotient [0, 1] /S„_i, which 
is a simplex and also has uniform measure. Then, taking the 
differences of successive points to obtain the probabilities pj 
is a linear isomorphism, which also preserves uniform mea- 
sure.) Then, if each child's birth is represented by a star which 
is also at a uniformly random position in [0,1], the probability 
of the jth birthday is exactly pj, the length of the jth interval. 

We note that the ordering of the stars and bars determines 
the number of children with each birthday. We claim that these 
multiset choices are all equally likely, as if the children had 
been bosons (with no distinguishing state other than the date 
of birth). This is made clear if we equivalently choose n — 
l+k points independently from / all at once, and then choose 
a random subset of n — 1 points to be the bars and the other 
k points to be the stars. These ((^)) = ( n ~l +k ) equally likely 
choices exactly correspond to a multiset choice of k unlabelled 
children distributed among « days, as claimed. □ 

We conclude with a version of the birthday paradox for 
fermions. 

Theorem 7 (Pauli). Given k fermions in any state on the exte- 
rior power A k (Jf?), there is no chance that any two have the 
same birthday. 

We leave the question of an anyonic birthday paradox, in- 
cluding non-abelian anyons, as a topic for future work. 
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